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Abstract 

Sensitive boundary condition dependence (BCD) is reported in a con- 
vectively unstable system with noise, where the amplitude of generated 
oscillatory dynamics in the downstream depends sensitively on the bound- 
ary value. This BCD is explained in terms of the manner in which the 
co-moving Lyapunov exponent (characterizing the convective instability) 
decreases from upstream to downstream. It is shown that a fractal BCD 
appears if the dynamics that represent the spatial change of the fixed 
point includes transient chaotic dynamics. By considering as an example 
a one-way-coupled map lattice, this theory for BCD is demonstrated. 

Convective instability, which causes amplification of a disturbance along 
flow|jl|, ^, is generally important in open flow systems ||. If a system is con- 
vectively unstable (CU), a small disturbance at an upstream position is amplified 
and transmitted downstream. Due to this property, spatiotemporal structure 
with a large amplitude can be generated in the downstream by a tiny fluctua- 
tion in the upstream. Such structure is referred to as noise-sustained structure 
(NSS)@. 

Convective instability is quantitatively characterized by a co-moving Lya- 
punov exponent A^, i.e., the Lyapunov exponent observed in an inertial system 
moving with the velocity w |^. If maxy Xy is positive for a given state, the 
state is convectively unstable. This condition is compared to that for linear in- 
stability, implies Aq > 0. Absolute stability (AS), which implies stability along 
any flowQ, ^, is guaranteed by the condition maxv Xy < 0. (With respect to an 
attractor, the co-moving Lyapunov exponent is used as an indicator of chaos: 
chaos with convective instability is characterized by the positivity of maXy Xy. 
However, with respect to a state, it is generally used to characterize its stability.) 

In a system with convective instability and noise, we have shown that the 
downstream dynamics depend on the boundary condition at the upstream. Such 



*fujimoto@complex.c.u-tokyo. ac.jp 



behavior exhibiting some threshold-type dependence on the boundary condition 
was identified and analyzed in connection with the change of the convective 
instability along the flowQ. In the present letter, we demonstrate that sensitive 
BCD of the downstream dynamics can appear in a class of noisy open-flow 
systems, and we clarify the condition for this appearance. 

In this letter we consider the simple case of a discrete system in one spatial 
dimension of the kind described above. As a simple example, we adopt a one- 
way coupled map lattice (OCML) §, |, |, |, [lo[ with noise: 

Here n is a discrete time step, i is the index denoting elements (i = 1, 2, , , , M = 
system size), and 77 is a white noise satisfying << 77^77^ >>= cr'^SnmSij, with 
<< ... >> representing an ensemble average (a <^ A fixed boundary 

condition is adopted and the (sensitive) dependence of the downstream dy- 
namics on a;° is studied. The use of a CML here is just for convenience for 
illustration. The results and theory we present are straightforwardly adapted 
to the case of coupled ordinary differential equations (ODE). 

In the present case we choose the logistic map f{x) = 1 — ax^ (—1 < x < 1). 
The parameters a and e are chosen so that in the noiseless case all elements 
are attracted to fixed points = x^ for any initial and boundary conditions 
[i.e., the fixed points (in time) are AS in the downstream]. This attraction to 
fixed points is realized in the strong coupling regime (see Ref.||^). The values 
of the fixed points can depend on the lattice site number (i.e. can be functions 
of space), and set of points {x^} form a spatially periodic pattern, as show in 
Fig.|l|. When noise is added, however, it can be amplified in the downstream to 
create oscillating motion (see Fig.||) if the upstream fixed points are convectively 
unstable ||]. 

Whether or not this noise-sustained structure is formed in the downstream 
depends on the boundary value xq and the noise strength. As a rough measure 
for the amplitude of the downstream oscillation, the root mean square (RMS) 
y'< a;^— < a;* >>^ is computed, where < ... > is the temporal average. As 
shown in Fig.^j, the RMS has a threshold- type dependence on xq- Such a BCD 
has been observed in a one-way coupled ODE||^. The mechanism responsible for 
the sensitive BCD clarified in that study is universal and can be summarized as 
follows: A change in x'~' causes a change in the value of upstream fixed points xl, 
which in turn causes a change in the degree of convective instability. Accordingly 
the downstream dynamics, generated through the spatial amplification of noise, 
can be different for different values of x'^ . 

In the present case this mechanism is quantitatively analyzed as follows: The 
spatial fixed points dynamics are described by the spatial recursive equation 

xl = (i-e)f{xl)+ef{xl-^) 



__ -l + \/l+4a(l-e){l-a<i(xi-^)^) (o) 
— 2a{l~e) ^ ' 

while the co-moving Lyapunov exponent Xv{i) of each fixed point xl ^ is 
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given by 

A.W = (l-e)log|/'(x:)|+dog|/'(xri)| 

+iog|i5fi+iogi^r- ^' 

A relevant quantity, characterizing the amphfication of a disturbance per lattice 
site is given by the spatial instability exponent X^{i) 

(i) = maXy (4) 

In Fig.|l|, A'^(i) is also plotted as a function of the lattice site number. For both 
the boundary values x° — 0.52 (□) and a:" = 0.60 (o), convective instability 
exists in the upstream, as indicated by X^{i). (Here the downstream pattern 
has spatial period 2, and X^{i) also oscillates with this period. For this reason, 
Sj=o'^^(* — ^)/2, the average over the spatial period, is also plotted.) For 

= 0.60, X^{i) (or the average over the spatial period) becomes negative at 
a smaller lattice site number than the case for x^ — 0.52. In fact, as discussed 
below, this difference in the convergence rate of X^ {i) is relevant to the BCD of 
the downstream dynamics. 

Since we have assumed that the fixed point pattern is AS in the downstream, 
the noise has to be amplified at a lattice point where the fixed point remains CU, 
in order for NSS to be formed. First, we estimate the lattice point defined as 
the site where the convective instability is lost; in other words, the site where the 
fixed point x\ changes from CU to AS. Recall that the x\ approach a periodic 
pattern in the downstream in the absence of noise. Denoting this spatial period 
by L, the lattice point is given by the point i such that X^fco^ X^ {i — i)/ L 
is positive for i < and negative for i > iu- As can be also expected by 
considering Figs.|l| and Fig.||, is strongly correlated with the relaxation scale 
of the spatial map. If the convergence of the spatial map to its attractor is more 
rapid, i„ is generally smaller. (See Fig.l: for example, iu = 12 for x° = 0.60 
indicated by •, while i„ — 36 for x^ = 0.52 indicated by x.) 

On the other hand, the scale ig required for the amplification of a tiny noise 
to 0(1) can be estimated by crexp(^^^j^ X^{i)) ^ 1. Then, the condition for 
the formation of NSS is simply given by 

iu-ig > 0. (5) 

In Fig.^, we have plotted iu — ig and the RMS of downstream dynamics as 
function of the boundary value x'^. The numerical results clearly support the 
conclusion that the condition for NSS is given by iu — ig > 0. Although this 
condition concerns only the sign of iu — ig, the amplitude of the NSS is also 
highly correlated with the value of iu — ig, as shown in Fig.^. In the case of 
Fig.^, the downstream AS dynamics are of spatial period 2, where NSS appears 
around = 0.52 and 0.97. The BCD here is simple, with just two regions 
allowing for NSS. For the case of spatial period 4, as shown in Figj^, there are 
many undulation, and this BCD has a self-similar fine structure, to be shown as 
fractals (see the blow-up of Fig.0) . The complexity of this self-similar structure 
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of the BCD increases with the period, as shown in Fig.^ for the case of spatial 
period 16. With this self-similar structure, a small difference in the boundary 
value results in a large difference in the downstream dynamics. 

We now discuss the origin of such sensitive BCD. As seen from Figs.^ and 
||, this BCD is due to the complicated structure of the BCD of z„ — ig- Here, 
ig has a rather smooth dependence on and the complicated structure is due 
mainly to z„. In fact, there are (infinitely) many local maxima of considered 
as a function of in analogy to the plot of z„ — ig in FigJJ. 

Recall that the scale i„ is highly correlated with the duration of the transient 
process of the spatial map Eq.(^), i.e., the number of steps required for an orbit, 
generated by the map starting from a:°, to fall into a periodic attractor. When 
the period of the attractor of the map g{x) is L, there are L stable fixed points 
for the map x g^{x). Each stable fixed point corresponds to a different 
phase of the periodic attractor of the spatial map g{x). For each value of the 
initial condition in the spatial map, the fixed point to which the g^{x) map 
is attracted [i.e., the phase of the cycle in the map g{x)\ is different. 

When multiple attractors coexist, it is often the case that the basin structure 
for each attractor is fractal iQ. This is true for g^{x), and here there can be 
infinitely many basin boundary points. For example, in Figj6| we have plotted 
the basin boundary points in the spatial map. For the spatial period 2 case, 
there are 4 such points, as shown in Fig.|^(a), while there are infinitely many 
points and the basin boundary is fractal for the case of spatial period 4 (or 
higher), as shown in Fig.||(b). Note that at unstable fixed points of g^(x), the 
phase of the attracted cycle slips. Successive preimages of such fixed points 
are nothing but the basin boundary of the map g^{x), with which fractal basin 
boundary is formed. 

For each basin boundary point, the number of transient steps before the 
attraction of an orbit to an attractor diverges. Hence, i„ takes a local max- 
imum at each point. Accordingly, there are infinitely many local maxima of 
i„, organized in a self-similar manner. Now, the sensitive dependence on x^ is 
understood resulting from a fractal basin boundary in the spatial map. Since 
this formation of a fractal basin boundary is rather common in one-dimensional 
maps (with topological chaos), sensitive BCD is expected to be a general phe- 
nomenon. Here, it should be stressed that although the mechanism here is based 
on the fractal basin in the spatial map, the sensitive dependence is explained as 
a dependence on the boundary condition rather than the initial condition! 15 . 

In conclusion, we have demonstrated the existence of a sensitive BCD in 
a system characterized by convective instability in the upstream. Note that 
the loss of convective instability in the downstream (for the noiseless case) is 
necessary to have such BCD. If the downstream without noise is CU, then the 
BCD becomes weaker as the distance from the boundary increases, and in an 
infinitely large system, this BCD eventually dies away completely. 

The origin of boundary condition sensitivity is found to be the complex 
transient dynamics in the spatial map, which lead to a fractal basin boundary. 
Accordingly, a sensitive BCD is generally expected as long as the spatial map 
exhibits topological (or transient) chaos. The analysis presented here can be 
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extended to the case in which, in the absence of noise, the downstream does not 
possess fixed points, but, rather possesses a stable cycle. 

Although we have studied the sensitive BCD for the simplest case with 
a CML, our analysis can be straightforwardly extended to systems of ODE. 
Convective instability in open flow is observed in a wide variety of systems with 
fluctuations, including chemical reaction networks optical networks 
traffic flow, and open fluid flow. Sensitive BCD is expected to be observed in an 
such systems. In particular, sensitive BCD in chemical reaction networks may 
be important in understanding diverse responses in signal transduction systems 
of cells. 

This work is partially supported by Grants-in-Aid for Scientific Research 
from the Ministry of Education, Science, and Culture of Japan (11CE2006 and 
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Figure 1: The patterns x\ and A'^(i). The upper figure displays the fixed- 
point pattern x\ in the absence of noise, while the lower figure displays the 
spatial pattern of Two sets of plots are overlaid for different values of 

the boundary, = 0.52 (solid line with □) and 0.60 (dotted line with o). For 
the lower figure Y^^^q A^(i — tj/2 (x for x^ = 0.52, • for 0.60) is also plotted. 
This difference reflects the difference in the downstream dynamics in the case 
with noise. For the boundary value for x!^ = 0.60 no oscillation is generated in 
the downstream (since the CU region is narrow, as discussed in the text), while 
oscillation is generated in the downstream for = 0.52, as shown in Fig.3. 
Here, a = 1.8, e = 0.83. 




Figure 2: Spatiotemporal plot of NSS with cr = 2 x lO'^, a = 1.8, e = 0.83, a;° = 
0.52. Here is plotted for every second lattice site. 
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Figure 3: BCD of the RMS for the downstream dynamics. NSS appears around 
x° = 0.52 and 0.97. Except these two regions for the boundary value, xl^ in the 
downstream falls on fixed points. Different lines correspond to different noise 
intensities: a = 2 x 10^^,2 x 10"'^ 2 x 10""*, 2 x IQ-^. With the increase in 
the strength of the noise, the size of the region in which NSS is found increases. 
Here e = 0.842, and the fixed point pattern is spatially period 2. 
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Figure 5: BCD of the RMS of x^^o with e = 0.928198, spatial period = 16, and 
a = 2 X 10~^. Successive blow-ups are given from the top to bottom figure. 
Each figure represents 5000 data points taken at equal intervals. 
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Figure 6: The orbits falling onto unstable fixed point for the spatial map g^{x). 
All the points on the orbits arc the basin boundary points, where the phase of 
the spatial cycle shows slips. There are 4 such points in the spatial period 2 
case (a), and an infinite number in the spatial period 4 case (b). 



